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We present updated results for the leptonic decay constants f K and fx, the light u, d, and s-quark 
masses, and the neutral kaon mixing parameter from mixed-action lattice simulations with 
staggered sea quarks and domain-wall valence quarks. We use the publicly-available 2+1 fla- 
vor MILC asqtad-improved staggered gauge configurations with multiple light sea-quark masses 
and three lattice spacings, and compute the kaon mixing matrix element with several partially- 
quenched valence-quark masses. We then extrapolate to the physical light-quark masses and the 
Q-i continuum using partially-quenched chiral perturbation theory formulated for mixed-action lat- 

r~ | tice simulations. For Bk we match the lattice four-fermion operator to the continuum using the 

nonperturbative method of Rome-Southampton. Our new results benefit from two significant im- 
t-H provements over our published work: (1) we have added a third lattice spacing of a«0.06 fm 

^ ! to better control the continuum extrapolation, and (2) we have implemented a new lattice renor- 

malization scheme (the RI/SMOIVL scheme developed by Sturm et al.) that suppresses chiral- 

00 

symmetry breaking and other infrared effects and, in practice, also shrinks the size of the 1-loop 
perturbative coefficient needed to match to the continuum MS scheme. When combined with the 
t-H use of volume-averaged momentum sources and twisted-boundary conditions, this significantly 

reduces the systematic uncertainty in the renormalization factor Zb k . 
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Figure 1: Global fit of the CKM unitarity triangle [3]. 

1. Motivation 

Lattice-QCD calculations of pseudoscalar decay constants, light-quark masses, and other kaon 
weak matrix elements are important ingredients for understanding the Standard Model and for 
constraining physics beyond the Standard Model. For example, the u, d, and s-quark masses are 
parametric inputs to calculations of Standard Model and new physics processes. The ratio of pseu- 
doscalar decay constants fa/ fit, when combined with experimental measurements of the leptonic 
decay rates, allows a precise determination of the ratio of CKM matrix elements |V„rf|/|Vu.s| [1]. 
The neutral kaon mixing parameter Bk, when combined with experimental measurements of indi- 
rect CP-violation in the kaon sector (Sk) constrains the apex of the CKM unitarity triangle. Finally, 
penguin-dominated K — > Tin and K — > tivv decays may be particularly sensitive probes of new 
physics once lattice weak matrix element calculations are sufficiently precise. 

Lattice QCD has played a key role in establishing that the CKM paradigm of CP-violation 
describes experimental observations at the ~10% level [2]. Many new-physics scenarios, however, 
predict new interactions between quarks and non-standard CP-violating phases. Given sufficient 
theoretical and experimental precision, these would lead to inconsistencies between measurements 
that are expected to be the same within the Standard Model CKM framework. Recent improve- 
ments in lattice weak matrix element calculations, especially of Bk and the fi-mixing SU (3)- 
breaking ratio ^, have revealed a ~3a tension that may indicate the presence of a non-Standard 
Model source of CP- violation (see Fig. 1) [4]. Given this tension, it is crucial to continue precision 
studies of kaon physics using multiple methods including our mixed-action approach. 

2. Overview of mixed-action calculation 

Table 1 shows the parameters of our numerical mixed-action lattice simulations. We analyze 
the 2+1 flavor asqtad-improved staggered gauge configurations generated by the MILC Collabora- 
tion [5]. Use of this large suite of ensembles with multiple spatial volumes ~ (2.5 - 4 fm) 3 , multiple 
light sea-quark masses ~ m s /\0 - m s , and three lattice spacings ~ 0.06 - 0.12 fm gives us good 
control over the systematic uncertainties associated with the chiral-continuum extrapolation. We 
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Table 1: Sea-quark ensembles and valence-quark masses used to obtain the preliminary results presented in 
this work. Ensembles shown in bold are new since our 2009 Bk publication [6]. 



generate domain-wall valence quarks at several partially-quenched masses ~ m s j 10 - m s . We apply 
HYP-smearing [7] to the valence domain-wall action in order to reduce the size of explicit chiral 
symmetry breaking [8] (m res ss 3 MeV at our coarsest lattice spacing and is ~ 30 times smaller 
at our finest lattice spacing). The use of domain-wall valence quarks makes the chiral-continuum 
extrapolation more continuum-like [9, 10] and simplifies the nonperturbative operator matching via 
the method of Rome-Southampton [11]. Hence our mixed-action approach is well-suited for weak 
matrix element calculations, as we show empirically in the next section. 

3. Testing the mixed-action method: decay constants and quark masses 

We compute the decay constant using a Ward identity to relate the axial current to the pseu- 
doscalar density; this method has the advantage that the renormalization factor is unity up to correc- 
tions of &(am KS ) ~ 10 3 — 10 5 . Once we have the decay constant for all combinations of valence- 
and sea-quark masses, we perform a combined chiral-continuum extrapolation via a simultaneous 
correlated fit of the full data set using NLO SU (3) mixed-action #PT [9] supplemented with higher- 
order analytic terms. Because the kaon mass is so heavy and our data is so precise, these terms are 
needed to interpolate about the strange-quark mass and obtain good confidence levels. We correct 
the numerical data for the known one-loop finite volume effects, and estimate the systematic un- 
certainty due to the chiral-continuum extrapolation by varying the fit function. As shown in Fig. 2, 
the value of f % in the continuum limit at the physical quark masses agrees with experiment. Our 
preliminary results for the pseudoscalar decay constants and their ratio are 

fn = 130.53(0.87) stat (1.68) zPT (0.80) F v(0.93) n (0.25) m9 MeV (3.1) 
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Figure 2: Chiral-continuum extrapolation of the pseudoscalar decay constant. Only data points with de- 
generate valence-quark masses are shown, but the fit includes all nondegenerate points. The circles denote 
"coarse" (a 0.12 fm) data, the squares denote "fine" (a 0.09 fm) data, and the triangles denote "superfine" 
data. The cyan band shows the continuum full QCD curve. The black "x " shows f % at the physical ud quark 
mass, where the inner solid error bar is statistical and the outer dashed error bar is the total systematic. For 
comparison, the red star shows the experimental result for f K , with which we are in good agreement. 




f K = 156.8(1.0) stat (l.l)^ PT (0.6) F v(0.8) ri (0.8) mf; MeV (3.2) 
hlfn = 1.202(0.011) stat (0.009) xF r(0.008)Fv(0.002) ri (0.005) w? , (3.3) 

where the errors are from statistics, the chiral-continuum extrapolation, finite volume effects, the 
scale uncertainty, and the light-quark mass uncertainties, respectively. We convert lattice quantities 
into physical units with the value of the scale r\ obtained by the HPQCD collaboration from several 
quantities including pseudoscalar decay constants and T splittings [12]. 

We determine the bare average ud and s domain-wall valence-quark masses by requiring 
that the pion and kaon masses in the continuum at these quark masses take their experimen- 
tal values [13]. We compute the quark-mass renormalization factor Z m = l/Z$ using a partly- 
nonperturbative approach in which we combine a nonperturbative determination of Za with a one- 
loop tadpole-improved lattice perturbation theory [14] calculation of (1 — Za/Zs). Because the 
ratio Za/Zs is close to unity, the size of the one-loop correction is small; nevertheless, the trunca- 
tion error is still the largest source of uncertainty in Z,„. We multiply the bare-quark masses by Z m 
to obtain the continuum MS masses; our preliminary results are 

94.2(1.4) stat (3.2) sys (4.7) raatch MeV (3.4) 
3.31(0.07) stat (0.20) sys (0.17) match MeV (3.5) 
4.73(0.09) stat (0.27) sys (0.24) match MeV (3.6) 



m MS(2 GeV) 
MS(2 GeV) 

m ud 

MS(2 GeV) 
m d 
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MS(2 GeV) = L9 o(o.08) stat (0.21) sys (0.10) match MeV (3.7) 
(m s /m ud ) = 28.4(0.5) st at(1.3) sys (m u /m d ) = 0.401 (0.013) stat (0.045) sys , (3.8) 



where the errors are from statistics, all systematics except for renormalization, and renormalization, 
respectively. Although we do not include isospin-breaking in our lattice simulations, we obtain m u 
and m c i separately following the method adopted by MILC [5], which uses the difference between 
the K + and K° meson masses and continuum estimates for the violation of Dashen's theorem. 

Our preliminary determinations of the decay constants and light-quark masses agree well with 
Nf = 2 + 1 lattice QCD calculations by other collaborations [15] (see also the recent review by 
Wittig [16]). This gives us confidence that the mixed-action method is reliable and can be used to 
obtain Bk and K — > %% matrix elements precisely with controlled systematic uncertainties. 



4. The neutral kaon mixing parameter: Bk 

In these proceedings we focus on improvements to our published calculation of Bk in Ref. [6]. 
Since 2009 we have addressed the two largest sources of uncertainty in our earlier work: 

1. The chiral -continuum extrapolation: We have added ensembles with larger volumes and 
lighter pions to shorten the chiral extrapolation. We have also added data at a third finer 
lattice spacing to reduce taste-breaking and shorten the continuum extrapolation. 

2. The AS = 2 operator renormalization: We now compute the renormalization factor Zb k us- 
ing volume-averaged momentum sources to reduce statistical errors and twisted boundary 
conditions to eliminate scatter from 0(4) symmetry breaking [17]. We also now compute 
Zb k with non-exceptional kinematics to reduce chiral symmetry breaking [18]. 

We compute Zb k via the nonperturbative method of Rome-Southampton [11]. We use an in- 
termediate lattice scheme with symmetric momentum p\ = p\ = (pi — p-i) 2 developed by Sturm 
et al. [18]; this suppresses chiral symmetry breaking and leads to a more convergent perturbative 
series for the conversion factor between the intermediate lattice RI/(S)MOM scheme and the con- 
tinuum MS scheme. Specifically, we use the RI/SMOM rti scheme (which has the same projectors 
as the RI/MOM scheme), for which the 1-loop correction to the continuum MS scheme at 2 GeV is 
~ 4 times smaller than for the standard RI/MOM scheme [19]. Because the RI/SMOM rti scheme 
is defined for massless quarks, we compute the bilinear and four-fermion vertex functions for sev- 
eral values of the valence, light, and strange sea-quark masses and extrapolate to the chiral limit 
tri^y \mi,m,h 0. The largest source of uncertainty in our calculation of Zb k is the residual pertur- 
bative truncation error from conversion to the continuum MS scheme. We estimate this in several 
ways including the size of 1-loop term (0.5%), comparison with lattice perturbation theory (2.0%), 
and the residual slope in {r\p) 2 (2.2%), and choose the largest error estimate to be conservative. 

We compute the Bk matrix element using symmetric and antisymmetric linear combinations 
of periodic and antiperiodic B.C. quark propagators to minimize finite temporal size effects. We 
obtain ~0.5-1.5% statistical errors in the lattice matrix element. We use the same approach for 
the chiral-continuum extrapolation as for the pseudoscalar masses and decay constants. We fit 
the renormalized Bk data to the NLO Si/ (3) mixed-action #PT expression [10] supplemented 
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Figure 3: Chiral-continuum extrapolation of Bk- The fit includes all available data, but only the nondegen- 
erate data points in which one valence-quark mass is close to the strange quark and the other valence-quark 
mass is the lightest on that ensemble are shown. The black " x " shows Bk in the continuum at the physical d 
and s quark masses, where the inner solid error bar is statistical and the outer dashed error bar includes the 
chiral extrapolation systematic. The star with dotted error bar, slightly offset, shows the total error in Bk- 



by higher-order analytic terms (see Fig. 3). We estimate the systematic uncertainty due to the 
chiral extrapolation by varying the fit function in several ways including using an analytic function 
without chiral logarithms, adding higher-order analytic and logarithm terms, and changing the 
value of f n in the coefficient of the 1-loop chiral logarithms. We conservatively take the largest 
difference from our preferred fit as the error. 

We obtain the following preliminary result for Bk in the MS scheme at 2 GeV: 

B^ S(2GeV) =0.5572(28) stat (45) xP T(33) F v(39) ri , mf/ (6) EM (134) match , (4.1) 

where the error labels are the same as in previous equations; the total uncertainty is 2.8%. The 
largest contribution to the error is from the renormalization factor; we aim to improve the uncer- 
tainty in Zb k due to the chiral extrapolation with additional data before publication. Ultimately, the 
dominant truncation error in Zb k may be reduced with a 2-loop continuum perturbative QCD calcu- 
lation. Our new result is ~ 1.1a higher than our published value, and agrees well with calculations 
from several other lattice methods [20] (see also the recent review by Mawhinney [21]). 

5. Summary and outlook 

Mixed- action lattice QCD simulations are well-suited to the calculation of weak matrix ele- 
ments, as shown by agreement with experiment and independent lattice results for pseudoscalar 
decay constants and light-quark masses, and also by our precise determination of Bk- Given this 
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and other recent improvements in lattice QCD calculations of Bk, it is no longer the dominant 
source of uncertainty in the Sk band. Hence, although recent Bk results slightly reduce the tension 
in the global CKM unitarity triangle fit, the tension remains at the 2-3 a level [22]. Given our 
success with Bk, we eventually aim to pursue K — > %% decays with the mixed-action approach. 
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